Constraining Light Colored Particles with Event Shapes 
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Using recently developed techniques for computing event shapes with Soft-Collinear Effective 
Theory, lep event shape data is used to derive strong model-independent bounds on new colored 
particles. In the effective field theory computation, colored particles contribute in loops not only to 
the running of a s but also to the running of hard, jet and soft functions. Moreover, the differential 
distribution in the effective theory explicitly probes many energy scales, so event shapes have strong 
sensitivity to new particle thresholds. Using thrust data from ALEPH and OPAL, colored adjoint 
fermions (such as a gluino) below 51.0 GeV are ruled out to 95% confidence. This is nearly an 
order-of-magnitude improvement over the previous model-independent bound of 6.3 GeV. 



Despite the fact that particle physics experiments have 
been running at and above 91 GeV center of mass ener- 
gies for over two decades, it is not known if the standard 
model represents the complete particle content below this 
scale. For particles which carry no standard model quan- 
tum numbers, the only hope of producing them at col- 
liders is through the Higgs, if there is a Higgs, and if 
they couple to it, or indirectly through off-shell inter- 
mediate states. But surprisingly, even colored particles, 
which interact with the strong force, are not significantly 
constrained. As long as they have small or vanishing cou- 
plings to electroweak gauge bosons, current data allows 
mass ranges well below the weak scale. A good example 
is a color adjoint Majorana fermion, such as the gluino 
in supersymmetric theories. 

The gluino is a color octet and thus should have a 
large production cross section at hadron colliders, a non- 
negligible contribution to four-jet events at lep, and a 
significant effect on the running of a s . Most of the cur- 
rent bounds on a color octet fermion depend on how it 
hadronizes and how it decays. For example, 

• If the gluino decays to two quarks and a very light 
neutralino, hadron collider data rules it out up to 
308 GeV at 95% confidence level (C.L.) Q. A re- 
cent study has shown that the Tevatron could probe 
gluino masses up to 150 GeV in the same decay 
channel independent of the neutralino mass 2] . 

• If the gluino is stable on detector lifetimes, ALEPH 
has excluded masses lighter than 26.9 GeV [3j. 

• A bound of 12 GeV, for a fixed a s (m z ) = 0.118, 
has been set based on the gluino's potential contri- 
bution to the parton distribution functions |^|. A 
strict lower bound (i.e., independent of a s ) has not 
been set. 

As for a model independent limit, ALEPH ,5] performed 
an analysis on four-jet observables as a measurement of 
the strong coupling constant and QCD color factors. The 
analysis found a good fit to QCD and ruled out gluinos 
below 6.3 GeV at 95% C.L.. An independent study [6j, 
which included the use of electroweak precision data, ar- 
rived at the same lower limit. In both cases, the bound 



comes essentially from the cross section for qqgg produc- 
tion which is very sensitive to the gluino mass at LEP 
1 energies. The scale 6.3 GeV is where these searches 
lose sensitivity and can be taken as the current model- 
independent bound on the gluino mass. 

Besides real production, new colored particles can be 
seen indirectly by their virtual effects. For example, 
any particle with color will contribute to the running of 
a s . Some of the current model- independent bounds come 
from fitting the 1-loop /^-function coefficient - which is 
sensitive to the number of flavors, n/ - to values of a s 
measured at different energies. For example, Delphi 
has done a study of mean values of event shapes and 
other inclusive observables leading to n/ = 4.7 ± 1.2 (and 
nf — 4.75 ± 0.44 when combined with low energy thrust 
data) and ruling out gluinos less than 5 GeV. This study 
includes data from 14 to 200 GeV. However, by averaging 
the observables - for example, into the mean thrust T(Q) 
- this approach is not optimized to take full advantage 
of the available data. 

A nominally positive feature of totally inclusive ob- 
servables, such as T(Q), is that only one scale appears, 
so the perturbation series in a s cannot be spoiled by the 
appearance of large logarithms. However, in searching for 
new particles through radiative corrections it is precisely 
these logarithms which have the most valuable informa- 
tion. In order to trust a differential calculation where the 
logarithms are relevant, the logs must be resummed. For 
many years resummation of event shapes was only avail- 
able at next-to- leading order [7J, which was insufficient 
to provide strong bounds on new physics because of large 
theoretical uncertainties. Recently, however, the thrust 
distribution was resummed to next-to-next-to-next-to- 
leading logarithmic order using techniques of effective 
field theory [Sj. Including matching to recent next-to- 
next-to-leading fixed-order (NNLO) event shapes [9], the 
theoretical uncertainty on the a s extraction from LEP 
was finally reduced to be sub-dominant to other uncer- 
tainties for the first time. Moreover, besides reducing the 
uncertainty, the effective theory approach makes explicit 
that a s is probed at many scales and so the sensitivity 
to new physics should be strong. Thus, it is natural to 



try to improve the model-independent bounds on new 
colored states using these recent theoretical advances. In 
this letter, we use these insights to improve the model- 
independent bound on the gluino by nearly an order of 
magnitude. 

The thrust distribution was shown in [a, [l0| to have 
the form 
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where r = 1 — T. Here, the matching function r(r) is de- 
fined as the difference between the fixed-order thrust dis- 
tribution and the fixed-order expansion of the resummcd 
distribution. We use r{r) at next-to-next-to-leading or- 
der, i.e. to a 3 s . 

The function R 2 (r) in Eq. ([1]) is the resummed distribu- 
tion. It can be calculated using Soft-Collinear Effective 
Theory Jill Il2l . Il3 | with insights from from [14J, |_15|, |lg, 
The result is 0: 

R 2 {t) = exp[4S(iJ. h , y) - 2A H (fi h , //) -8S(fi ,u) 
+AA J (n j ,n) + 4S(n s ,n) + 2A 3 (ji.,n)] 
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The derivatives in Eq. @ are to be taken analytically and 
then n set to its canonical value n = 4Ar(/Xj, Ms)- Here, 
H(Q, (/,) is the hard function and j(£, /i) and s(£, m) are 
the Laplace transforms of the jet and soft functions; all of 
these have power series expansions. S(u, m) and Ax{v 1 m) 
are auxiliary functions defined as integrals over various 
anomalous dimensions. Explicit expressions for these 
functions can be found in [8(] . The scale jjl in Eq. ^ is ar- 
bitrary and the distribution is formally independent of it, 
but different values of /1 can be chosen for calculational 
convenience. 

The formula ([2]) is a simplified version of the one in [8| , 
valid when the scales are set to their canonical values: 



Hh = Q, Mj = QVt, Ms = Q T 



(3) 



As mentioned above, a calculation of mean thrust, or a 
fixed-order calculation of differential thrust, would only 
probe a s (Q) at a single scale, the hard scale \xu = Q- 
But the differential thrust distribution probes even lower 
scales. For example, in the two-jet limit thrust reduces 
to the sum of hemisphere masses, Q 2 t ~ M\\ + M^. 
The effective theory expression associates this mass scale 
with the scale of jet functions, and probes it through 
Mi ~ Q\/t. Actually, the effective theory makes it ap- 
parent that even lower scales, associated with soft modes 
of QCD, are relevant. These are probed by the soft scale 
Ms ~ I 1 "]/ Hh ~ Qt, which is a type of seesaw scale lower 
than both of the physical external scales Q and Qy/r [10(. 
Since a s is larger at lower energy, resumming logs of the 



soft scale is critical to generating an accurate thrust dis- 
tribution. 

Because the differential thrust distribution is sensitive 
to many scales, it would be sensitive to the presence of 
new colored particles with a variety of masses. These new 
states would affect the running of a s , through the QCD 
beta function, as well as the hard, jet, soft anomalous 
dimensions - which appear implicitly in ^ - and the 
fixed-order hard, jet, and soft functions, H,j and s". 

Throughout the following we modify the standard 
model by adding Any new flavors of mass m at a thresh- 
old scale Mth- For example, a new massive quark cor- 
responds to Anf — 1 and a gluino to Anj = 3 [20|. 
Below the scale Mth> the new flavors are integrated out, 
insuring decoupling as m — > 00. This will, in general, in- 
duce discontinuities in a s (n) and in the hard, jet and soft 
functions, all of which are unphysical by themselves. The 
resulting thrust distribution, however, must be smooth. 
In fact, one can show that the effective field theory dis- 
tribution is independent of pt t h order-by-order in pertur- 
bation theory 2l|. For simplicity, we take Mth = ti and 
match a s at one-loop. To avoid having to run the jet and 
soft functions through the threshold, we choose fi = m in 
Eq.([2]) when fi s < m < fij. For m < fj, s , we take fi = fi s 
and for m > fj,j we take [i = [ij. 

To demonstrate the sensitivity to new states, we be- 
gin by looking at a single data set, the ALEPH data from 
lep 1 at 91.2 GeV [22|. We perform a bin- by-bin correc- 
tion for hadronization and quark masses using PYTHIA 
v.6.409 [ll]. Using the fit region 0.10 < r < 0.24 the 
soft scale fi s probes 9-22 GeV and the jet scale Mj probes 
29 - 45 GeV. Thus, if there are n f flavors below 9 GeV, 
we do not have to worry about an explicit threshold and 
may simply run a s and the other objects using this value 
of n/ throughout (a more refined procedure is described 
below). To derive a bound on the number of light flavors, 
we perform a least-squares fit to the experimental data. 
For the errors used in the fit, we include both the ex- 
perimental statistical uncertainty and also the statistical 
uncertainty in the fixed-order thrust distribution. The 
fixed-order result was calculated numerically, with some- 
what slow convergence at NNLO, and to be conservative 
we rescale the NNLO uncertainties by a factor of 1.5 to 
account for the fact that the errors may have been un- 
derestimated. A combined fit with two free parameters 
gives a s (m z ) = 0.1169 ± 0.0004 and n f = 5.32 ± 0.59, 
where the errors are statistical only. 

For a second example, using the same LEP 1 data set, 
we note that a gluino of mass m = 25 GeV lies outside the 
range of scales probed by the hard, jet, and soft functions. 
Thus, it can be modeled by taking Anf — 3 for the hard 
and jet functions, and Anf — for the soft function. 
Performing the fit with these values, we find x 2 = 31.7 
with the gluino compared to x 2 = H-9 for the standard 
model, with 13 degrees of freedom. The fits for the two 
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FIG. 1: Theoretical prediction versus ALEPH data at lep 1 
for the standard model and the standard model with a 25 
GeV gluino. The total statistical uncertainty band includes 
theoretical statistical uncertainty from the Monte Carlo used 
to generate the NNLO fixed-order thrust distribution. 
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FIG. 2: Bounds on light colored particles from lep data. The 
darker region is completely excluded at 95% confidence. The 
lighter region is an uncertainty band including estimates of 
various theoretical uncertainties. 



models are shown in Figure [TJ where it is clear that the 
model with the gluino is systematically worse. 

To properly scan over masses, we must specify how 
to handle the thresholds. First, consider the total 
hadronic cross section, Chad- The exact leading or- 
der dependence of Chad on the new particle mass can 
be extracted from [28]. For m < /i, the contribution 
to the total cross section is proportional to A<7 n ad = 
o2(A«) (M$) + Pr($) + \ log(^)) , where p v is the 
virtual contribution which vanishes at m = oo and pn 
is the real emission contribution which vanishes for m > 
Q/2. The explicit log compensates the //-dependence of 
a s and is necessary to have a smooth m — > limit. We 
will use this exact expression Acrhad for the new physics 
contribution to Chad in Eq. (JTJ), but observe that, as 
shown in [28], it is well approximated for < m < Q 
by the leading power in m 2 /Q 2 . Actually, it is not 
clear whether the experiments would have included de- 
cay products of real gluinos in their event selection for the 
thrust distribution, so in the spirit of providing a model- 
independent bound, we allow A<7had to scan between 
and the cross section for An/ additional massless fla- 
vors. This variation is included in the uncertainty band 
described below. 

The exact contributions of massive colored states to 
the jet, soft, and hard functions are not known, but since 
the same loops and real-emission diagrams are relevant 
for them as for Arjhadj it is likely that the result would 
be similar to that of Aa^^. Thus, we assume the leading 
power is linear in m 2 / p? h for the hard function, m 2 / \ih for 
the jet function, and m 2 / p 2 for the soft function. That is, 
we take H, j and s" to interpolate between the expression 
for rif = 5 + An/ flavors at m — and n/ = 5 flavors 
at the relevant threshold. This removes any remaining 
discontinuity in the thrust distribution, and should be a 



good approximation to the (unknown) exact result. In 
a similar vein, the matching correction, r(r) in Eq. ^, 
formally takes place at the hard scale Q. However, it 
depends on n/ and would be discontinuous as m crosses 
Q unless the discontinuity is removed by inclusion of ex- 
plicit mass corrections. We use an interpolation also lin- 
ear in m 2 /Q 2 for this effect. Using this model for the 
mass thresholds, in lieu of the exact result, introduces 
some theoretical uncertainty. To account for that un- 
certainty, we explore some variations of the model and 
include the errors in our final bound, as described below. 

With this treatment of the threshold effects, the thrust 
distribution is smooth and can be compared with the 
data for each m and An/ . We perform a combined fit to 
the ALEPHJp and opal [3, [13] data sets from 91.2-206 
GeV [26|, |237 The fit regions used are 0.1 < r < 0.24 
for lep 1 , and 0.04 < r < 0.25 for ALEPH lep 2 and 
0.05 < t < 0.22 for opal lep 2 . The data are cor- 
rected bin-by-bin for hadronization and bottom/ charm 
mass effects using PYTHIA. We perform a least-squares 
fit of the theoretical prediction to the corrected data, us- 
ing errors which include both the experimental statistical 
errors and the statistical errors of the NNLO fixed-order 
calculation, rescaled by 1.5, as described above. For the 
standard model, the x 2 is 85.7 for 78 degrees-of-freedom. 
For each value of m and An/, we minimize x 2 an d com- 
pute the maximum likelihood ratio as compared with the 
standard model. The resulting 95% C.L. bound is shown 
in Figure [2] For An/ = 3, the limit is ra-g > 52.5 GeV. 
For a real gluino (with the appropriate group theory fac- 
tors differing from An/ = 3 at higher orders), the bound 
differs by 0.03 GeV. 

To account for the theoretical uncertainty, we include 
an uncertainty band (the light shaded region in Figure[2|). 
This subsumes the following variations: (?) Removing 
the lowest bins from each data set in the fit. (it) Not 



interpolating the total cross-section and matching cor- 
rection (we include variations both with rif = 5 and 
with nt = 5 + Any in <7i ia( j and r(r)). (in) Varying pa- 
rameters in the hadronization model between PYTHiA's 
default values and the ALEPH [29j and opal [30] opti- 
mized tunings (P AR J (81) = 0.290, 0.292, 0.250 and 
PARJ(82)=1.00, 1.57, 1.90). (iv) Using a power cor- 
rection proportional to m/[i x instead of m 2 / ' \j? x for the 
threshold corrections. The band in Figure [2] includes the 
maximal and minimal bounds at 95% C.L. for each value 
of to and Any. For the gluino, this gives m-g > 51.0 — 54.0 
GeV. For our final bound, we take the least restrictive 
value, nig > 51.0 GeV. 

From Figure [2j it is clear that this method has the 
strongest sensitivity in an intermediate mass range, 10 
GeV < to < 40 GeV. This range roughly coincides with 
the scales probed by the jet and soft functions in the 
fit regions of the thrust distributions. For masses below 
about 10 GeV, the mass threshold lies outside of the fit 
regions and the effect on the event shape can be partially 
compensated by a change in a s . When the mass falls 
inside the range of the thrust distribution, it is more dif- 
ficult to compensate by rescaling a s , hence the stronger 
bound. With additional independent constraints on a s , 
for example, from the lattice [3l| or from r decays [32], 
one might be able to close the light mass window more 
tightly. This might, for example, even rule out additional 
light colored triplets or scalar adjoints. However, as the 
lattice and r-decay determinations of a s (which take 
place at similar scales) are themselves inconsistent by 
more than two standard deviations, it is unclear whether 
a definitive bound could be obtained in this way. The 
main result of this paper is that event shapes alone are 
sufficient to exclude light and intermediate mass gluinos 
up to 51 GeV, independently of their decays. 
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